day21

mo = AOI::aoi_get("Mount Saint Helens") %>% 
  AOI::aoi_buffer(.5)

elev  = elevatr::get_elev_raster(mo, z = 12)

writeRaster(elev, "/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif", overwrite = TRUE)

wbt_slope("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")

wbt_aspect("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_aspect.tif")

wbt_hillshade("/Users/xingxin/Github/geog176a-summer-2020-lab1/mount_sthelens-elev.tif","/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_hillshade.tif")

elevt = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")

slop = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_slope.tif")
  
aspe = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_aspect.tif")  

hisd = raster("/Users/xingxin/Github/geog176a-summer-2020-lab1/mst_hillshade.tif")

plot(elevt, axe = F, box = F, col = viridis::viridis(256), main = "elevation")

plot(slop, axe = F, box = F, col = terrain.colors(256), main = "slope")

plot(aspe, axe = F, box = F, col = rainbow(8), main = "aspect")

plot(hisd, axe = F, box = F, col = gray.colors(256, alpha = .8), main = "hillshade")